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Abstract 



We present an exact, unconstrained representation of the electron opera- 
tors in terms of operators of opposite statistics. We propose a path-integral 
representation for the t-J model and introduce a parameter controlling the 
semiclassical behaviour. We extend the functional approach to the Hubbard 
model and show that the mean-field theory is equivalent to considering, at 
Hamiltonian level, the Falikov-Kimball model. Connections with a bond- 
charge model are also discussed. 
71.20.Ad, 75.10.Lp, 74.65.+n. 



Typeset using REVT^ 



The investigation of the nature of the ground-state and of the low-lying excitations of 
strongly correlated electron systems is a fundamental problem of the modern many-body 
theory. The standard model of correlated electrons on a lattice is the Hubbard Hamiltonian 

Hh = -tJ2 ^^j4,sCj,s + (1) 
i,j,s i 

where Ajj is the symmetric adherence matrix connecting the nearest neighbor (n.n.) sites 
of a hypercubic lattice of M sites, cf^ and Ci^s are creation and annihilation operators for 
electrons with spin projection s =t, J., respectively, and ni^s = cJ,sCi,s. Henceforth we denote 
the states at site i with the notation 

|l,^)i = 4^|z/), |0,^), = |z/)i, 

|1,4). = 4Jz/), |0,4). = 44k), (2) 

where \v) = |z/)j. The essential behaviour of Hh in the strong-coupling regime U/\t\ ^ 1 
below half-filling is commonly investigated [|I| by means of the effective Hamiltonian H'^^ = 
Htj + H^^^ acting in the subspace without doubly occupied states |0,JJ.)j, where 

Hu = -tY: AAsC,,s + Jj:iS. S, -\n,N,), (3) 

is the t-J model and H^^\ not explicitly displayed |Q, is the three-site term. Here 
denotes the summation over the n.n. sites, J = At^/U, Ci^s = Ci,s(l — ni^g), with s denoting 
the opposite spin projection, are projected electron operators, and 5^ = 1 Y.s,p '^\.s^spCi,p and 
c\ Jii^s are the spin and number operators, respectively, with cr^p the Pauli matrices. 
To approach the problem of strong correlation, slave-boson (-fermion) methods |l|,|^ have 
been widely employed, because they give a physically intuitive way to work in the subspace 
without doubly occupied states and allow the introduction of mean-field approximations by 
assuming condensation of the bosons. In the present investigation we provide a path-integral 
description of the models (0) and (^, starting from the observation that the decomposition 
of the electron operators via operators of opposite statistics can be achieved in terms of 
exact operator identities. Hence, contrary to the slave-particle methods, in our approach 
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constraints among the operators are absent. In fact, the states (||) can be generated by 
means of canonical spinless fermions f- and Pauh operators aj 

11,^). = ^.,-!^), |o,^), = /M,-l^^), (4) 

where crj^± = ((7j_^±icri^y)/2 and |f) = <S>i 1^)1 is the reference vacuum annihilated by the /j's 
and CTj^+'s. We observe that the representation is closed, since no new states can be generated 
by applying fj and/or crj^_, so that the operators of a given basis can be expressed in terms 
of the others 

where Ui^s = 1 — "^2,5 and 7i,± = (1 ± crj^2)/2. Hence, it is possible to construct nonlinear 
combinations of the electron operators satisfying to {//, fj} = 6ij, [crj_a, dj^is] = 2i 5ije"-^^aj^^ 
or, relevant for us, one sees that the anticommutation relations {c| ,,, Cj^p} = Sij6sp can be 
fulfilled via operators of opposite statistics. 

For the moment, it is useful to avoid referance to a particular Fock representation. 
To this aim we say that the model (|T]) [or (^] acts on "objects" |^, cr)j of four (or 
three) different species, placed exactly one to a site, which we label by defining the "grade" 
7r(^) = (— )^ = ±1 = 0, 1), i.e., by dividing the species into even and odd, and the 
"isospin" = ±1 (pictorially cr =lt, J|) quantum numbers. The unitary transformations of 
the species assigned to each site are given by Hubbard operators Xf^,^, = cr)i{^f, af\i, 
which form the fundamental representation of the su(2|2) graded algebra The special 
character is due to the completeness relation J2^a ^i,Co- = 1, where Ni^^„ = Xf^^ are the local 
densities, expressing the one-to-one correspondence between objects and sites. With our 
notation, the anticommuting, or "odd" operators, are those changing the grade [i.e., with 
7r(,^) = — 7r(^/)], henceforth denoted 

— If^'^ i/f^ — Vlf^ T- — y'^^ — (7o\ 
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so that Xi = and = rl^. For the commuting, or "even" operators, [i.e., those with 
7r(^) = 7r(^/)] we define the hnear combinations 





1 ^ Y'^^ 

2 2 




- 




Si,t = 


2 i'^t 2 




Y^ii- 

— 


(7b) 


Li^z = 


2"^iofr 2 








Li^t = 


1 Y^t 1 -'- 

2"^iofr 2 






(7c) 



For graded (i.e., supersymmetric) algebras, the grading of the states is a convention of 
purely formal nature, because the grading of the operators does not depend on the choice 
used. Physically, this freedom corresponds to the fact that the particle-hole transformation 
|0,cr)i ^ leaves su(2|2) unaltered, i.e., 

{Si^t, Si) <-> (Li t, Li), 
ixt, Xi, 1"/, ^ {Xi,t, Xi,i^, n,i}.), (8) 

where Si^± = Si^^ =t iSi^y and Lj = Lj^^. ± iLi^y. Hereafter we shall always use the corre- 
spondence (1^), so that Si and Lj are identified as the local spin and pseudospin operators, 
respectively. Henceforth we also define the isospin vector as the sum of spin and pseudospin: 
Qi = Si + Li. The total numbers of species-(^, a) objects N^a = J2i ^j,Co- are related to the 
total numbers of spin-s electrons, henceforth denoted A^^, through the self-evident relations 
N-^ = Ni-f^ + A^ojj. and A^'j^ = A^'ijj, + A'qjj., from which we obtain the equivalent expressions 

Sz = liN^^~N,^^), Sz = liN^-N^), 

L, = ^{No^-No^), Lz = ^{M-N^~N^), (9) 

where S = Y^i^i and L = Y^i^i- The operators Sz and L^ are conserved because N-^ and 
are constant, however from Eq. (^) we might also say that conservation follows because 
even and odd objects, when not conserved, are both created and/or destroyed in pairs. The 
latter is a useful way to visualize the less intuitive properties of the pseudospin. The picture 
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is sharpened by observing that spin and pseudospin always act on A^^o- as raising-lowering 
operators, thus changing by one unit in states with a definite number of objects. We 
have [Nq.^, L±] = ±L±, and [A^o^, L±] = tL±, as well as the commutators where Nq^j Ni^ 
and L± ^ S±. 

Hubbard operators are useful for discussing symmetry properties and by using them it is 
easily seen that any transformation of the electron operators can be exactly rephrased in the 
basis (0). In particular, in the representations (|]) and (H) the particle-hole transformation 
(I) is achieved by letting c|_| ^ and // ^ fi, respectively. The realization of Eq. (7) in 
the basis (|) is 

x! = li,+fi, = (^i,+fi, Si^t = ^(1 - rii), 

Xi=(^i,-fi, rl=li,-fi, Si = ^ai{l-ni), (10) 

where Ui = f} fi, and the operators not explicitly displayed are easily obtained by means 
of Eq. (g). The local densities are Ni^^ = 'ji^+rii, Ni^^^ = -fi-Ui, Ni^i^ = 7i,+(l - n^), 
~ 7«,-(l ~ "^i)) automatically satisfy to Z^go-^i.^o- = 1- For the realization of Eq. 
(7) with electron operators we refer instead to Ref. 

We now consider the auxiliary model acting by permutation of the four species of 
objects |^,cr)i 



(11) 

{id) ' 

where Plj^ = 2{SiSj + Si^tSj^t) and = 2{LiLj + Li^Lj^) are the spin and pseudospin 
permutators acting nontrivially on the odd and even objects, respectively, and P^^ = XiXj,a+ 
T^Tj^„ + h.c. is the operator permuting pairs of objects of opposite grade. Eq. (|llD, henceforth 
referred to as the extended model, is the case x = y = 2 oi the class of Hamiltonian iy(^l^) 
permuting objects of y odd and x even different species, introduced by Sutherland for the 
special case J = ±2|t| in studying exactly solvable systems in one dimension. The spectrum 
of if contains that of //(^l^) by construction and noting that the latter is actually the t-J 
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model, we see that information about the Hamiltonian (^) can be obtained from the model 
(|n|) by imposing the conserved constraint Nq^^ = 0. Using Eq. (^) iJ^^I^) becomes 

HS = t E A.,P.,/;/, + ^ E - 1)' (12) 

where Pij = |(1 + (Tj (X,) and Ajj = {1 ~ rii — rij), and where we have added a constant to 
ensure that Eq. (|T2]) reduces to the standard definition (|) of the t-J model for A'o^ = 0. 
We consider the generalization (p!T| ) useful because it leads in H^J both to the isotropic 
contribution Pij in the hopping term, similarly to the generalization of KhaliuUin 0], and 
a) to the quadratic contribution in the magnetic term. In this respect, Eq. (0) is easier 
to study than the original t-J model, as for, e.g., the magnetic term of Hfj in the basis @ 
has a four-fermion interaction Mij = (1 — ni){l — rij). 

For the time being the state |0,JJ.)j will be interpreted as a fictitious "polarization state 
of the hole" , and not as the physical doubly occupied state, because we are interested in Eq. 
(|l^) only as an auxiliary description of Eq. (j^). From this interpretation it follows that the 
conserved quantity G = Y^a ^la always corresponds to the number of electrons N^i = J2s 
and 6 = J^a^Oa/M = to the doping. The extended model commutes with S and L, 

and from Eq. (|^) one easily sees that for any given G = Nd the condition A^o^ = projects 
onto the sector where pseudospin and total pseudospin L attain their maximum allowed 
value L, = L = {M - Nd)/2. 

In a grand-canonical approach the partition function of the model ([I^ ) is 

Z = TTexp{-pH^}, H^ = H^^-J2fi"^No^, (13) 

where are the chemical potentials for the two species of "holes" . In the basis a path- 
integral representation for Z can be built-up immediately by using Grassmann variables 
rji, 7]* for the spinless fermions and standard SU(2) coherent states for the isospin vectors. 
In fact, Qi = Si + Li is the sum of two reducible operators, for which vector addition does 
not apply. The only eigenvalues of Qi^a are = ±^ and indeed in the basis (H) we can write 
the isospin vector as a pure spin-| operator: Qi = \(7i. Moreover, this property allows us to 
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introduce an expansion parameter, denoted g, by enlarging the dimensionality of the SU(2) 
representation of Qi. For consistency we enlarge the whole even sector (|7^,c) of su(2|2), 
setting Si^t = ?(1 — ni) and Li^t = Qni in Eq. (p!0|). Following the spin-wave approach 
adopted in Ref. 0, the generalization of at arbitrary q is thus achieved by letting 

(P^^._l)^p{-)=2(Q4--g'), (14) 
P,, ^ P/+) = 2(4Q, + ^^^^ ^ q ± Q^^^. 



Due to the presence of in Eq. (|T2|), the Grassmann integration over the variables rj* 
and rji entering Z can be performed exactly and we obtain Z = J VQ exp{—Ses}, where VQ 
denotes the usual integration over classical unit-vectors i7j(r) and the effective lattice 
action is 

^efi = ^qT.f drA^dA + 4g'^ E f dr vlj^ (15) 



-Tr ln{[9. - + /i_fi.,.)]5., + VtA,,p{+^} 

with Ai the Dirac potential for a monopole of unit strength, V-^^ = ^{QiQj ± 1), and 

The model (|T^) behaves as a correlated hopping, and in fact /ij enters Ses as a rotationally 
invariant local chemical potential correlated to the background. In dealing with the extended 
model H^J itself, the most natural choice is to treat the two species of holes equally, thus 
setting /i_ = 0, whereas when considering the t-J model one has to impose the condition 
din Z /dfi^^ = 0. As expected, in the latter case, the rotational symmetry in isospin space is 
explicitly broken, as for /i_ is coupled to Qi^z and acts as an effective magnetic field. However, 
as discussed in Ref. 0, in the static limit t = the ground-state energies of H^J and Htj 
are degenerate and we expect the degeneracy to remain for a physically sensible range of 
the parameters J/\t\ and 6. Hence, in this case, one can investigate the zero temperature 
properties of the t-J model directly with the extended Hamiltonian H^j, thus working with 
the simpler condition /x_ = 0. 



The doping is accounted for by the third term in the r.h.s. of Eq. (pISl). This contribution 
is ineffective at zero doping, and because for 6 = spin and isospin coincide (as for Lj = 0), 
then ^eflf correctly reduces to the standard lattice action of a quantum antiferromagnet [H. 



We also notice that, contrary to other proposed functional integrals |jTO[ for the t-J model, 
the imaginary part of Eq. (|T5|) is a sum of standard Berry phases, so that the quantization 
condition on the monopole flux ^ is satisfied for any 6. We remark that the classical field 
i7j(r) is the expectation value of the operator Qi/q, hence for 6^0 physical correlation 
functions can be evaluated only by differentiating Z with respect to suitable source terms, 
e.g., Hs = J2ihiSi = J2ihi{l — ni)Qi, added to H^. Allowing for these modifications of 
Eq. (|13D, the effective action is modified as well. However, the resulting Eq. (^) remains a 



purely classical expression, for which a gradient-expansion [|ri|] controlled by the parameter 
q allows extraction of the relevant low-energy behaviour. 

The basis (H) is suited to define an effective action also for the Hubbard model. To 
this aim one only needs to insert Eq. (|^) into Eq. ([1|). However, before performing this 
straightforward calculation, we divide the lattice into two sublattices A and B, and perform 
the transformation c|,, — > e*^'c|,,, where 9i = for i E A and Oi = ^ for i G -B, so 
that Hh is mapped into an equivalent Hamiltonian Hh, whose hopping matrix satisfies to 
tAij = —tAji = tA*j. We make this change because it is known ||l2l that the operator M, 



with components M_ = 4a + M+ = m1, and = - c[iCi,i)' where 

rrii = (1) for i G A (B), commutes with H = Hh + ULz- By changing the phases we 
have M Q, because the alternating sign = (— )™'' cancels, and Hh Hh, so that 
[Hh + ULz,Q] = 0. Hence, this choice of the phases leads to the representation of the 
Hubbard model 

HH = -tY. h [P^lflf3 + {P^1 " 1)^^] + ^at, (17) 

where = §(/,/,■ - /]//) and H^t = UE j7j_?T,i, with a manifestly isospin invariant 
hopping. In Eq. (|T^ one easily identifies a term Hq conserving the number of objects of 
each species, and a term Hj leading instead to transitions AG = ±2. Setting Aij —Aij via 



the transformation // — > e*^*//, we have Hq ^ Hq = Hs,t + tJ2i,j-^ijPijfifj^ and confronting 
this expression with the results in Ref. [§] one sees that Hq is also the particular bond-charge 
model i^bc = + tJ2i.j,s-^iM,s(^j,s{^i,s + ^j,s)- The equivalence of these different looking 



models is found because both have the same expression Hq = i/at — i Y.(i,j) Pij in terms 
of abstract Hubbard operators, and holds because a global phase change in one basis leads 
instead to site-dependent transformations when rephrased in the other, e.g., 

fl^e'^fl ^ 4,^expM2n,,,-l)}4. (18) 

This is actually the main reason why, before making use of Eq. (^, we have brought Hh to 
a form with antisymmetric adherence matrix. By replacing directly Eq. into Hh and 
M, one obtains messy representations obscuring the result [Hh + UL^, M] = 0. 

Using Eq. (|I4|) and following the discussion leading to Eq. (0), the Hamiltonian H^ = 
Hh — P-Nei can be extended at arbitrary q and its partition function represented as the 
path-integral Z = J VQexp{—S^}, 

S^s = + f drA^OA - ^ In det [iT.f ] , 

where /xq = yu/(2g) and, setting Uq = U/{2q), h = fiQ — Uo/2, 

Kll = [dr + Aq\^ + hn,,;)] 5,, - AqHk,,v\l\ 

Klf = -AqHK,v\;^ , KJl = -AqHK,v\;^ , (19) 



r|f = [dr - V(^ + /^a,.)] 5., - ^qHk,v^\ 



The square root of the determinant of the block-matrix K^^ enters into the effective action 
because [Q, due to the presence of Cij = ^{ViVj ~ VjVi) the classical counterpart of Hh, 



we have performed the integration with respect to the variables r/f, with r]} = rji, rjf = rj*, 
satisfying to the condition {r]^)* = y]i(Jxfia- The Pfaffian structure of S'Jf is quite natural: 
The off-diagonal blocks vanish only for a ferromagnetic background, and setting VLi^a = ^az 
one can easily see that in this case both effective actions Ses and S*^ correctly reduce to 
the logarithm of a free spinless fermion determinant. At half-filling one has = 0, so that 
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is rotationally invariant in the isospin. The symmetry is instead exphcitly broken away 
from half-fiUing, as expected because of the relation qfli^z = (^il'S'j^^ + — Nci^i)\^i) , where 
iVei^i is the local electron number operator. 

In general temporal fluctuations of the classical fields fij(r) (i.e., quantum fluctuations of 
the Qi's) are suppressed by letting g oo. Assuming the generalization (0), in Eq. one 
has c| | oc q±Qi^z and cj^ oc Qi-, and at large isospin we consistently [l^ obtain c- ^ = 0{q) 
and = O(^), so that the propagation of the spin-| electrons becomes smaller and 
eventually vanishes at g = oo. Hence, at mean-field level reduces to the effective action 
S t hat one would have obtained by considering the Falikov-Kimball Hamiltonian [jl],|l6| 



This is a useful result, because one can reasonably obtain sensible information about the 
Hubbard model by systematically including corrections starting from a model, the Falikov- 



Kimball, for which many exact results are known [16 
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